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Abstract 

Background: Biofortification of staple crops is a cost effective and sustainable approach that can help combat 
vitamin A and other micronutrient deficiencies in developing countries. PCR -based DNA markers distinguishing 
alleles of three key genes of maize endosperm carotenoid biosynthesis {PSY1, IcyE and crtRBl) have been developed 
to facilitate maize provitamin A biofortification via marker assisted selection. Previous studies of these functional 
DNA markers revealed inconsistent effects. The germplasm previously employed for discovering and 
validating these functional markers was mainly of temperate origin containing low frequencies of the favourable 
allele of the most significant polymorphism, crtRB1-S f JE. Here, we investigate the vitamin A biofortification potential 
of these DNA markers in a germplasm panel of diverse tropical yellow maize inbred lines, with mixed genetic 
backgrounds of temperate and tropical germplasm to identify the most effective diagnostic markers for 
vitamin A biofortification. 

Results: The functional DNA markers crtRB1-5'TE and crtRB1-3 r T\i were consistently and strongly associated with 
provitamin A content across the tropical maize inbred lines tested. The alleles detected by these two functional 
markers were in high linkage disequilibrium (R 2 = 0.75) and occurred in relatively high frequency (18%). Genotypes 
combining the favourable alleles at the two loci (N = 20) displayed a 3.22 fold average increase in (3-carotene 
content compared to those genotypes lacking the favourable alleles (N = 106). The PSY1 markers were monomorphic 
across all of the inbred lines. The functional DNA markers for IcyE were associated with lutein, and with the ratio of 
carotenoids in the alpha and beta branches, but not with provitamin A levels. However, the combined effects of the 
two genes were stronger than their individual effects on all carotenoids. 

Conclusions: Tropical maize inbred lines harbouring the favourable alleles of the crtRBI -S f JE and 3'TE functional 
markers produce higher levels of provitamin A. Such maize lines can be used as donor parents to speed up the 
development of provitamin A biofortified tropical maize varieties adapted to growing conditions and consumer 
preferences, providing a route towards mitigation of vitamin A malnutrition in Sub-Saharan Africa. 
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Background 

Carotenoids are naturally occurring organic pigments pro- 
duced by plants mainly as integral component of the light 
capturing and protective plastidal apparatus [1]. Carot- 
enoid intake plays an important role in human nutrition 
and health owing to the association of their consump- 
tion levels with reduced risk of diseases such as cardio- 
vascular disease, cancer, and age-related sight problems 
arising from deficiencies of lutein and zeaxantine [2,3]. 
Provitamin A carotenoids including a-carotene, p-carotene 
and p-cryptozanthine are precursors of vitamin A, a micro- 
nutrient essential for normal development and functioning 
of the human body [4,5] . Vitamin A deficiency is currently 
a global public health problem inflicting morbidity, stunted 
growth, night blindness and loss of both sight and lives 
in the developing world [6,7]. Vitamin A deficiency is 
estimated to affect 190 million preschool children and 
19 million pregnant and/or lactating women world wide 
[8], which aggravates poverty and underdevelopment chal- 
lenges in developing countries [9,10]. 

Maize can naturally accumulate both provitamin A and 
non-provitamin A carotenoids in its kernel, and is known 
for its genetic diversity of carotenoid content and profiles 
[11-13]. However, provitamin A usually constitutes only 
10 to 20% of the total carotenoids in maize kernel, and the 
commonly cultivated and consumed yellow maize cultivars 
have less than 2 \ig g 1 DW provitamin A [14]. Exploitation 
of the natural genetic diversity of maize in carotenoids 
through biofortification by combining conventional and 
molecular breeding can increase provitamin A concentra- 
tion in maize endosperm [15,16]. Such increases in pro- 
vitamin A concentration can be beneficial for public health 
in Africa where maize is a major security and staple crop 
for more than 300 million people [17,18]. Biofortification 
offers a safe, effective, cheap and sustainable approach to 
combating vitamin A and other micronutrient deficiencies 
[14,19-23]. 

One of the major challenges in maize breeding for 
high provitamin A levels is the quantification of caroten- 
oids in the endosperm of large number of breeding lines. 
High performance liquid chromatography (HPLC) is the 
commonly used method for carotenoids analysis be- 
cause of its accuracy. However, HPLC is expensive, time 
consuming and relatively low throughput, limiting its 
use for routine breeding within resource-limited plant 
breeding programs [24]. Even though the variation in 
intensity of yellow color in maize endosperm is attribut- 
able to the variations in carotenoid content and profile, 
selection for high provitamin A maize based on kernel 
color is not reliable due to its poor correlations with 
provitamin A content [12,25]. Marker assisted selection 
[26] using functional DNA markers [27] offers an effect- 
ive tool for screening a large number of breeding mate- 
rials for their carotenoid profile and content accurately 



and cheaply within a reasonable timeframe within a plant 
breeding program. 

The carotenoid biosynthesis pathway is well studied in 
plants [28-31], where the genes encoding the enzymes of 
the biosynthesis pathway are known [1,32]. Specific nu- 
cleotide sequence variants within the key carotenogenic 
genes have also been characterized, and shown to contrib- 
ute significantly to accumulation of provitamin A and total 
carotenoids in maize endosperm [12,33,34]. For instance, 
[12] showed that four polymorphic sites in the gene encod- 
ing lycopene epsilon cyclase (IcyE) were associated with 
the variation in ratio of carotenoids in the a to (3 branches 
of the carotenoid biosynthesis pathway, leading to a three- 
fold increase in provitamin A. [34] also identified three 
polymorphic sites in another downstream gene encoding a 
p-carotene hydroxylase enzyme (crtRBl) accounting for 
40% of the observed variation in p-carotene concentration 
in maize endosperm. Two polymorphisms in the gene 
encoding phytoene synthase (PSY1) have been identified 
explaining 7 to 8% of the variation in total carotenoids 
[33]. These findings allowed the development of breeder 
friendly PCR based functional DNA markers that can be 
used as tools for detecting alleles representing each of the 
polymorphic sites in the three genes. 

While these functional DNA markers can be used to 
facilitate the development of maize cultivars fortified with 
high provitamin A, their efficacy in breeding lines used for 
maize variety delivery to Sub-Saharan Africa has not been 
fully elucidated. Some studies have examined the individual 
and combined effects of the functional polymorphisms of 
IcyE and crtRBl on carotenoids [11,33,35,36]. It has been 
observed that the proposed diagnostic polymorphisms for 
IcyE could not distinguish between inbred lines with high 
lutein and high zeaxanthine, representing carotenoids in 
the a- and p-branches of the pathway, respectively [11]. In- 
bred lines having high p-carotene contents have been iden- 
tified, although they were carrying the unfavourable alleles 
of crtRBl -5 'TE and -3'TE. Similar inconsistencies for the 
favourable allele of the crtRBl -3TE marker have been ob- 
served [36]. A recent study tested two of the three signifi- 
cant polymorphic sites of IcyE (5TE and 3' indel) and one 
of the three functional polymorphisms of crtRBl (3TE) 
using 26 tropical segregating populations [35]. Their results 
showed that the effects of IcyE on both ratio of a to p 
branch carotenoids and total provitamin A content were 
inconsistent across the populations, whereas the crtRBl -3 ' 
TE polymorphic site had a large effect on p-carotene and 
provitamin A concentrations. In contrast, significant effects 
have been detected for all the functional polymorphisms 
for individual and haplotypes of selected polymorphisms 
of IcyE, crtRBl and PSY1, using inbred lines with trop- 
ical, subtropical and temperate backgrounds [33]. The 
reported inconsistencies in the effects of the diagnostic 
DNA markers for the proposed favourable alleles of 
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IcyE and crtRBl require further investigation of these 
markers using diverse inbred lines to identify the most 
robust and effective markers for marker assisted 
selection. 

In the present study, the functional DNA markers for 
IcyE, crtRBl and PSY1 were tested on a set of diverse 
tropical inbred lines with mixed genetic backgrounds of 
tropical and temperate origin developed within the 
International Institute for Tropical Agriculture (IITA) 
maize breeding program for Africa. The inbred lines were 
first evaluated for carotenoid content and composition 
across two seasons (in 2010 and 2011) and exhibited con- 
trasting variations. Two of the functional markers of crtRBl 
(i.e. the 5TE and 3TE markers) were found to be in high 
linkage disequilibrium, and displayed consistent and strong 
effect on the provitamin A carotenoid contents of the 
inbred lines. The deployment of these functional DNA 
markers can accelerate the biofortification of tropical- 
adapted maize varieties with elevated levels of provitamin A. 

Methods 

Plant materials 

One hundred and thirty diverse tropical adapted yellow 
maize inbred lines were assayed for carotenoid profiles 
and content and used in this marker-trait association study. 
These inbred lines were developed within the maize breed- 
ing program for Africa at the International Institute for 
Tropical Agriculture (IITA) from eight bi-parental crosses 
of tropical inbred lines, four broad based populations, and 
28 backcrosses involving temperate lines as donors of high 
p-carotene (Table 1). 

Field evaluation 

The 130 inbred maize lines were field evaluated at IITAs 
research site (7°29 '11.99 "N, 3°54'2.88"E, altitude 190 m) 
in Ibadan, Nigeria, in 2010 and 2011. The field trial was 
arranged in a 13 x 10 alpha-lattice design with two repli- 
cations. Each inbred line was planted in a single 5 m long 
row with spacing of 0.75 m between rows and 0.25 m 
between plants within a row. Different fields were used in 
each season. Fertilizer and field management practices 
recommended for optimum maize production were used. 
Seed samples for carotenoid analyses were produced by 
self pollination of at least 5 representative plants in each 
row. Self-pollinated ears in each row were harvested, dried 
under ambient temperature, and threshed, with minimal 
exposure to direct sunlight. One hundred kernels were 
drawn from seed samples for carotenoid analysis. 

Carotenoid analysis 

Carotenoids were extracted from the maize kernels and 
quantified by HPLC at the University of Wisconsin, USA. 
The extraction protocol and carotenoid analysis used was 
the method described in [37]. Briefly, 0.5 g finely ground 



sample of each entry was transferred into a 50 ml glass 
centrifuge tube to which 6 ml of Ethanol plus 0.1% butyl- 
ated hydroxyl toluene were added, vortexed for 15 sec- 
onds, and incubated in 85°C water bath for 5 min. 500 (il 
of 80% potassium hydroxide (w/v) was added to each sam- 
ple, vortexed for 15 seconds, and incubated in the 85°C 
water bath for 10 min with vortexing at about 5 min inter- 
val. Samples were then immediately placed on ice and 
3 ml ice cold deionized water added to each of them, vor- 
texed for 15 seconds, and 200 \A internal standard p- Ape- 
s' -carotenal and 4 ml hexane added. After vortexing and 
centrifugation, the top hexane layer formed was trans- 
ferred into a new test tube. The hexane extraction was re- 
peated twice, adding 3 ml hexane each time. Samples were 
allowed to dry down completely under nitrogen gas using 
a Turbovap LV concentrator (Caliper Life Sciences) and 
reconstituted in 500 \A of 50:50 MethanobDichloroethane. 

Fifty micro-liter aliquots of each extract were injected 
into an HPLC system (Water Corporation, Milford, MA). 
The Waters HPLC components was operated with Em- 
power 1 software and included a 717 Plus auto sampler 
with temperature control set at 5°C, Waters 1525 binary 
HPLC pump, and a 2996 photodiode array detector for ca- 
rotenoid quantification. Carotenoids were separated by C30 
Carotenoid Column (4.6 x 250 mm; 3 (im) eluted by a mo- 
bile phase gradient from 100% methanol/water (92:8 v/v) 
with 10 mM ammonium acetate to 50% methyl tertiary 
butyl ether. The flow rate was 1.0 mL/min and the solvents 
were HPLC grade. To maximize detection of carotenoids, 
absorbance was measured at 450 nm. Alpha-carotene, p- 
carotene (cis and trans isomers), p-cryptoxanthin, lutein, 
and zeaxanthin were quantified. 

Total carotenoid was calculated as the sum of concen- 
trations of a-carotene, lutein, p-carotene, p-cryptoxanthine 
and zeaxanthine. Provitamin A was calculated as the sum 
of p-carotene and half of each of p-cryptoxanthin and 
a-carotene concentrations, since the latter two contribute 
50% of the value of p-carotene as provitamin A [38]. Other 
derived carotenoid traits were also calculated as indicated 
in [12] and [34], namely the ratio of carotenoids in p to a 
branch of the carotenoid pathway, the ratio of p-carotene 
to p-cryptoxanthine and the ratio of p-carotene to total ca- 
rotenoids. The natural logarithms of the ratios were calcu- 
lated to allow statistical analysis of the data, as the ratios 
followed a highly non-normal distribution. All concentra- 
tions were described in \ig g" 1 dry weight (DW). 

PCR based genotyping 

For PCR based genotyping of the functional DNA markers, 
leaf samples were collected from 3 to 4 randomly selected 
plants of each inbred line within one of the replications of 
the field trial described above at 40 days after planting. 
DNA samples were isolated from freeze dried leaf sam- 
ples of each genotype using either a CTAB (cetyl trimethyl 
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Table 1 Maize genotypes used in present study 



Serial no. Origin Number 

of lines 

1 4001 1 

2 9450 1 

3 KU1409 1 

4 9450/KI28 2 

5 9450/KI21 6 

6 9450/CM 116/9450 2 

7 9450/KI21 -1 -4-1 -1 -1 -B/DE3/9450/KI21 -1 -4-1 -1 -1 -B 4 

8 9450/KI28-1 -2-1 -1 -B/DE3/9450/KI28-1 -2-1 -1 -B 3 

9 9450/KI21 -1 -4-1 -1 -1 -B/DE3/9450/KI21 -1 -4-1 -1 -1 -B 1 

1 0 9450/KI21 -1 -5-3-2-1 -B/DE3/9450/KI21 -1 -5-3-2-1 -B 1 

1 1 9450/KI21 -1 -5-3-2-2-B/DE3/9450/KI21 -1 -5-3-2-2-B 2 

1 2 9450/KI2 1 -3-2-2- 1 -3/KU 1 409/M0 1 7LP A/KU 1 409 1 

13 DE3/KU1414-SR/KU1414-SR 3 

14 KU1409/NC358/KU1409 3 

15 KU1414-SR 1 

16 KU1414-SR/KVI11 2 

17 KU1409/DE3/KU1414-SR 1 

18 KU1409/DE3/KU1409 23 

19 KU1409/KU1414-SR/A619 9 

20 KU 1 409/KU 1 41 4-SR/KV1 1 1 1 

21 KU1 409/KU 1 41 4-SR/KVI3 12 

22 KU1 409/KU 1 41 4-SR/M162W 3 

23 KU1 409/KU 1 41 4-SR/NC298 8 

24 KU1 409/KU 1 41 4-SR/NC350 10 

25 KU1 409/KU 1 41 4-SR/SC55 2 

26 KU1409/SC55/KU1409 4 

27 KU 1 41 4-SR/CI7/KU 1 41 4-SR 2 

28 KU1414-SR/CML328/KU1414-SR 1 

29 POP66SR/ACR91 SUWAN1 -SRC1/ACR91 SUWAN1- 2 
SRC 1 - 1 /SYN-Y-STR-34- 1 -1-1-1-2-1 -B*3 

30 POP66SR/ACR91 SUWAN1 -SRC1/ACR91 SUWAN1- 2 
SRC 1 -4/400 1 /KI2 1 -4- 1 - 1 -1 -1 

31 POP66SR/ACR91 SUWAN1 -SRC1/ACR91 SUWAN1- 3 
SRC1 -6/(MP420/4001 /MP420) 

32 POP66SR/ACR91 SUWAN1 -SRC1/ACR91 SUWAN1- 2 
SRC1 -8/POP61 -SR-1 1 -2-3-3-1 -B 

33 POP66SR/ACR91 SUWAN1 -SRC1/ACR91 SUWAN1- 1 



SRC1 -9/(9450/CM1 1 6/9450)-3-3-1 -2-1 

34 SYN-Y-STR-34-1 -1 -1 -1 -2-1 -B*3/(DE3/CI7)/SYN-Y-STR- 1 
34-1 -1-1-1 -2-1 -B 

35 SYN-Y-STR-34-1 -1 -1 -1 -2-1 -B*5/NC354/SYN-Y-STR-34- 1 



1 -1-1-1 -2-1 -B*5 

36 SYN-Y-STR 1 

37 4205/CI7/4205 1 

38 1 



Table 1 Maize genotypes used in present study 

(Continued) 

ACR9^TZL-CCOMP1 -Y-S3-1 3-1 -B*2/CI7/ACR9^ZL- 
CCOMP1-Y-S3-13-1-B*2 

39 ACR97TZL-CCOMP1-Y 2 

40 SC55/KU1414-SR/KU1414-SR 1 

41 TZE-COMP5-Y-C7 1 

42 Z.DipIo BC4 1 



ammonium bromide) based DNA extraction protocol or 
QIAGEN DNeasy® Plant Mini kit (Qiagen Inc., Hilden, 
Germany) following the company's protocol 

PCR based functional markers of three genes IcyE, 
crtRBl and PSY1 were deployed across all the 130 in- 
bred lines. PCR conditions, cycling profiles and primers 
used were based on those reported by [12] for IcyE [34], 
for crtRBl and [33] for PSY1. The primers used to amplify 
the lcyE-3 TE indel marker were forward 5-ACCCGT 
ACGTCGTTCATCTC-3 ' and reverse 5-ACCCTGCGT 
GGTCTCAAC-3' [35]. Primer oligos were obtained from 
Integrated DNA Technology Inc (IDT, Belgium). All PCRs 
were run using BIOTA Q™ DNA polymerase kit (Bioline 
Ltd, UK) with a mixture composed of 2 \A lOx NH4 PCR 
buffer, 1 \A of each primer, 1 or 1.5 \A (depending on the 
marker) of 50 mM MgCl 2 , 0.15 [A of BIOTAQ™ polymer- 
ase, 1 \A of Dimethyl Sulfoxide (DMSO) to enhance speci- 
ficity, and ultra pure water making up to 25 \A total 
volume. PCR fragments were confirmed by sequencing 
three samples representing each allele of the 6 functional 
markers. PCR product sequences were aligned with se- 
quences of the three genes, downloaded from GenBank 
of NCBI or MaizeGDB, using CLC genomics workbench 
(CLC Bio, Denmark) sequence analysis software. Frag- 
ments in the PCR products were resolved using 2% w/v 
super fine resolution (SFR™) agarose gel. The names of 
polymorphic sites of each gene and the nature of poly- 
morphisms are indicated in Table 2 according to their 
respective references. 

Statistical analysis 

The carotenoid data was analyzed using PROC MIXED 
procedure of SAS version 9.3 (SAS Institute, Cary NC) 
based on alpha lattice design in which lines were treated 
as fixed effects, while blocks, replications, years, and year 
by line interaction were treated as random effects. Esti- 
mates of repeatability were calculated as indicated in 
[33]. Letters for mean separation were generated using a 
SAS macro [39]. Spearman rank correlation coefficient 
was calculated using SAS 9.3 (SAS Institute, Cary NC) 
to test the consistency of ranking of the inbred lines for 
accumulation of carotenoids across seasons [13]. 

Associations between variation in carotenoid concentra- 
tion and markers of each gene were calculated using the 
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Table 2 Nomenclature of functional DNA markers and their allelic series [12,33,34] 



Gene 


Polymorphic site/marker 
gene name-polymorphism) 


Nature of polymorphism 


Allelic series and notations* 


PSY1 [33] 


PSY-SNP7 


A-C substitution SNP 


A,C 




PSY1-IDI 


378 bp indel 


0, 378 


LCYE[U] 


LCYE-5'TE 


285 indel 


1, 2, 3, 4 




LCYE-SNP (216) 


G-C SNP 


G,T 




LCYE-3 f Me\ 


8 bp indel 


8, 0 


crtRBl [34] 


crtRB1-5'TE 


397/206 bp indel 


1,2, 3 




crtRB1-lnDel4 


12 bp indel 


12, 0 




crtRB1-3'TE 


325/1250 bp indel 


1, 2,3 



^Allelic variants denoted in bold face underlined letters represent the best favourable alleles as proposed in the references. In the current study IcyE 5'TE yielded 
no amplification for 73% of the inbred lines invariably. Hence there is an additional notation in the results section for those samples scored as a '0' allele to mean 
'no amplification'. 



mixed linear model (MLM) [40] implemented in TASSEL 
version 3.0 [41]. MLM incorporates population structure 
and kinship in the analysis to control spurious association 
results [40]. Best linear unbiased estimates (BLUEs) calcu- 
lated via the generalized linear model (GLM) option by 
selecting only the phenotype data was used for association 
analysis combined across the two seasons [41], Linkage 
disequilibrium between functional markers was also 
calculated using the same software. Population structure 
(principal component analysis, PCA) and kinship of the 
130 inbred lines were estimated within TASSEL 3.0 
using 62,000 SNPs that covered the 10 maize chromo- 
somes generated by genotyping by sequencing (GBS) 
method at the Institute for Genomic Diversity (IGD), 
Cornell, USA, according to [42]. The SNPs were filtered 
out from the GBS pipeline output using a threshold of 
5% minimum allele frequency and 20% maximum missing 
data. In addition, 2328 SNPs were filtered using 0% missing 
and 5% minimum minor allele frequency and used for 
hierarchical clustering of SNP data for 26 inbred lines that 
harbored the favourable alleles of the two most significantly 
associated markers (crtRBl-5'TE and -3'TE) to assess their 



genetic diversity. The unweighted pair-group method with 
arithmetical averages (UPGMA) provided in PowerMarker 
version 3.25 [43] was employed to construct a dendrogram 
from Neis 1972 frequency based distance matrix [44], Sin- 
gle and two-way ANOVA were conducted to determine 
genotype effects of the functional polymorphisms using 
PROC GLM and PROC MIXED of SAS 9.3 (SAS Institute, 
Cary NC). 

Results 

Carotenoid profiles and levels are diverse across the 
maize inbred lines 

The analysis of variance (ANOVA) combined over the 
two years revealed highly significant variation among the 
maize inbred lines for all carotenoids, except for a- 
carotene (Table 3). The effects of year, and year by line, 
interaction were significant on all carotenoids, except (3- 
cryptoxanthine. High repeatability estimates ranging from 
62 to 89%, were recorded for all carotenoids expect for 
a- carotene demonstrating the importance of the genetic 
component of the total variation observed for the traits. 
Replication did not have a significant effect on all 



Table 3 Combined ANOVA for carotenoid content of 130 inbred lines evaluated in 2010 and 2011 

Sources of variation DF Mean squares of carotenoids 3 



lut zeax (3cry acar Pear pva tear 



Line 


129 


26.81*** 


133.99*** 


12.37*** 


0.43 


40.71*** 


4] 49*** 


226.72*** 


Year 


1 


519.32*** 


878.82** 


7.95 


51.98*** 


48.63** 


84.07** 


332* 


Rep (Year) 


2 


2.81 


22.07 


0.69 


0 


0.57 


1.29 


23.58 


Block (Rep*Year) 


36 


3.75** 


13.43*** 


0.44 


0.01* 


0.88 


1.36 


31.32*** 


Line*Year 


129 


10.32*** 


14.21*** 


1 .62*** 


0.37*** 


5.21*** 


6.63*** 


43.45*** 


Residual 


222 


2.15 


4.08 


0.48 


0 


0.76 


1.01 


11.28 



r 0.62 0.89 0.87 0.15 0.87 0.84 0.81 



a Carotenoids are abbreviated as lut Lutein, zeax Zeaxanthine, [Scry p-cryptoxanthine, acar alpha-carotene, (Scar p-carotene, pva provitamin A, tear Total carotenoid, 

r repeatability, DF degrees of freedom. 

* **, *** = significant at P < 0.05, 0.01, and 0.001, respectively. 
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carotenoids. Spearman's rank correlation coefficients 
across years were significant (p < 0.001) for p-carotene (r = 
0.83), p-cryptoxanthine, (r = 0.75) zeaxanthine (r = 0.86) 
a-carotene (r = 0.30) and lutein (r = 0.49). Zeaxanthine was 
the dominant carotenoid identified with average mean 
value of 9.66 \ig g" 1 followed by p-carotene, 4.21 \ig g 1 and 
lutein, 3.58 \ig g" 1 . The a-carotene contents of most of the 
maize inbred lines were very low and not significantly dif- 
ferent from zero (apart from 16 inbred lines). Estimated 
means averaged over the two years varied from 0.45 to 
13.51 \ig g' 1 for lutein, from 0.04 to 25.90 \ig g 1 for zeax- 
anthine, from 0.08 to 8.55 \ig g 1 for p-cryptoxanthine, 0 to 
16.38 \ig g" 1 for p-carotene, from 0 to 17.25 \ig g l for pro- 
vitamin A, and from 4.43 to 42.71 \ig g' 1 for total caroten- 
oids (Table 4, Figure la and b). 

Several of IITA's tropical adapted inbred lines harbour 
alleles of IcyE and crtRBl markers proposed for elevated 
provitamin A in maize endosperm 

The two markers of PSY1 [33] were monomorphic for 
the favourable allelic variants across all 130 inbred lines, 



and were thus not considered for further analysis. In 
contrast, all the PCR markers for IcyE and crtRBl were 
polymorphic across the inbred lines. Sequences of all sam- 
pled PCR fragments were aligned to their corresponding 
gene sequences (data not presented) confirming their iden- 
tity. Alleles 2 and 4 of the 5TE polymorphic site of IcyE, 
and allele 2 of the 3TE and allele 3 of the 5TE polymor- 
phisms of crtRBlwere not detected in this study (Table 5). 
Frequencies of the IcyE favourable alleles varied from 12% 
to 83%, while those of crtRBl varied from 18% to 19% 
(Table 5). 

The favourable alleles identified by the crtRBl -5'TE 
and crtRBl -3 'TE markers were present in 26 inbred lines, 
co-occurring in 20 inbred lines. The two polymorphisms 
showed high linkage disequilibrium (R 2 = 0.76). However, 
both (i.e. crtRBl -5'TE and crtRBl -3'TE markers) were not 
in linkage disequilibrium with the crtRBl -indel4 marker 
(Table 6). Linkage disequilibrium values between markers 
of IcyE and crtRBl genes were low with R 2 values ranging 
from 0.004 to 0.188. 

Five different donor lines were determined to have con- 
tributed the favourable alleles of crtRBl-5'TE and 3'TE to 



Table 4 Carotenoid levels and IcyE and crtRBl genotypes of selected maize inbred lines 



Carotenoids* (ng g _1 clry weight) Genotype^ 



Line entry no. 


lut 


zeax 


(kry 


acar 


|3car 


tpva 


tear 




IcyE 






crtRBl 




















5TE 


SNP(216) 


3' InDel 


5'TE 


InDel 4 


3'TE 


39 


170ab 


2.72c 


0.44c 


1.14a 


16.38a 


17.25a 


22.88abc 


'0' 


G 


8 


2 


0 


I 


45 


1.03b 


2.39c 


0.34c 


0.97a 


14.79ab 


15.45ab 


19.53bc 


'0' 


G 


8 


2 


0 


113 


124 


0.99b 


1.12c 


0.16c 


1.18a 


13.64abc 


14.33ab 


1 7.08bc 


3 


T 


0 


2 


0|12 


113 


106 


1 .80ab 


0.26c 


1.30c 


1.21a 


11. 80a be 


13.02abc 


16.30bc 


'0' 


G 


8 


2 


12 




99 


5.80ab 


2.32c 


0.39c 


1.20a 


11. 30a be 


12.1 labc 


21. 10a be 


3 


T 


0 


2 


0 


1 


23 


5.76ab 


4.83bc 


0.75c 


1.02a 


10.77abc 


11.63abc 


23.04abc 


'0' 


G 


8 


2 


12 


I 


101 


12.09a 


5.02bc 


0.83c 


1.39a 


9.76abcd 


10.91abc 


29.1 9abc 


'0' 


G 


8 


2 


12 


I 


107 


1 .90ab 


1.05c 


0.33c 


0.80a 


9.07abcd 


9.65abc 


13.17 c 


'0' 


G 


8|0 


2 


12 


I 


11 


2.03ab 


10.80abc 


1 .99bc 


0.83a 


7.61 bed 


9.06abc 


23.48abc 


'0' 


G 


8 




0 


3 


50 


7.28ab 


18.74a 


8.55a 


0.57a 


7.54bcd 


12.10abc 


42.71a 


3 


T 


8|0 




0 


3 


92 


3.96ab 


15.56ab 


6.02ab 


0.70a 


6.72 cd 


lO.IOabc 


32.98abc 


'0' 


G 


8 




12 


3 


98 


3.22ab 


9.43abc 


2.59bc 


0.77a 


6.30 cd 


7.98bc 


22.27abc 




G 


0 




0 


113 


120 


1.37ab 


10.36abc 


1 .92bc 


0.60a 


6.27 cd 


7.52bc 


20.46bc 


'0' 


G 


8 




0 


3 


96 


8.89ab 


19.82a 


3.28bc 


0.37a 


2.93d 


4.76c 


35.41 ab 


'0' 


G 


0 




0 


3 


SE 


1.95 


2.36 


0.67 


0.45 


1.23 


1.40 


3.53 














Min. 


0.45 


0.04 


0.08 


0 


0.03 


0.06 


4.43 














Max. 


13.51 


25.90 


8.55 


1.68 


16.38 


17.25 


42.71 














Grand mean 


3.58 


9.66 


2.92 


0.4 


4.21 


5.87 


20.78 















*Means within column followed by same letters are not significantly different using least square means and Tukey's multiple comparison test (p < 0.05). 
Abbreviations of carotenoids are described in legend for Table 3. 

**Allelic series of genotypes composed of six polymorphisms: IcyE 5'TE, IcyE SNP (216), IcyE 3'indel, crtRBl 5'TE, crtRBl lnDel4, and crtRBl 3'TE from left to right 
(discussed in later sections). Heterozygosity of the T allele of the dominant marker SNP (216) is not determined. Heterozygous alleles are separated by '\. 
Favourable alleles are bolded and underlined; description of the markers is presented in Table 2. The haplotype 3, T, 8|- 1, 0, 3 represents the worst allelic 
combination. The haplotype '0', G, 8|0, 2, 12, 1 represents the best allelic combination. 
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Figure 1 Distribution of mean concentrations of carotenoids for 130 inbred lines, a) Histogram; error bars represent standard error of least 
square means of the respective carotenoid concentration, b) Box plots; endpoints of upper and lower whiskers represent maximum and 
minimum concentrations, respectively; upper and lower edges of boxes represent third and first quartiles, respectively; line inside box represent 
median; symbol ♦, represent mean. Abbreviations of carotenoids described under Table 3. 



the 26 inbred lines carrying either one or both of these 
favourable alleles. The largest proportion of the inbred lines 
having favourable alleles of crtRBl-S'TE and/or -3 'TE were 
derived from backcrosses containing a temperate inbred 
line DE3 as a donor parent. These inbred lines are among 
those lines that exhibited the highest levels of provitamin 
A. In addition, three tropical lines were the recurrent par- 
ents of the best inbred lines, also carrying favourable alleles. 
Cluster analysis based on Neis 1972 frequency based dis- 
tance using UPGMA separated the 26 best "favourable- 
allele -carrying" inbred lines into three major groups with 
one line separated from the major groups (Figure 2). Even 
though lines originating from the same backcross were 
grouped together; they showed considerable levels of within 
group diversity. 



crtRBI functional markers had the largest effect on 
provitamin A variation across the 130 inbred lines 

Association analyses were conducted to determine the 
relationship between polymorphic alleles (for IcyE and 
crtRBI) and phenotypes (carotenoid levels and profiles) 
based on each year mean and means averaged over two 
years. Alpha-carotene was excluded as a phenotype from 
the analysis due to its extremely low concentrations and 
lack of significant variability among the maize inbred lines. 
The results of the association analysis are presented in 
Table 7 and Figure 3. The 3TE and 5TE polymorphic 
sites of the crtRBI candidate gene were found to be sig- 
nificantly associated with carotenoids, and all the derived 
traits, consistently over the two years (a = 0.01). The ex- 
ceptions were lutein levels which were not affected by 
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Table 5 Observed alleles and frequencies of the 
favourable allelic class of PSY1, IcyE and crtRBl functional 
markers 



Marker 


Expected 
allelic 


Allelic 
variants 

UUjCI vcu 


Favourable 
allele 


Frequency of 
the favourable 
allele (%) 


PSY1 SNP7 


A, C 


A 


A 


100 


PSY1 InDell 


0, 378 


0 


0 


100 


IcyE 5'TE* 


1,2, 3,4 


'0', 1, 3 


1 


12 


IcyE SNP 
(216) 


G, T 


G, T 


G 


83 


IcyE 3'indel 


8, 0 


8, 0 


0 


38 


crtRBl 5'TE 


1, 2, 3 


1, 2 


2 


18 


crtRBl indel4 


12, 0 


12, 0 


12 


19 


crtRBl 3'TE 


1, 2, 3 


1,3 


1 


18 



*For IcyE 5TE the vast majority of the inbred lines (> 70%) did not yield 
any amplification, and thus scored as '0' alleles to represent no 
amplification, not deletion. 

**There were individuals that were heterozygous for some of the marker 
classes. Description of markers and expected alleles are presented 
in Table 2. 



both crtRBl markers, and the a to |3 branch carotenoids 
ratio which was not affected by crtRB 1-5 TE in the second 
year. The two crtRBl markers explained from 13 to 53% 
of the variation in carotenoids and derived traits in the 
first year, and from 17 to 63% in the second year inferring 
from the R 2 values (Table 7). CrtRBl -indel4 accounted for 
9% of the variation in provitamin A (a = 0.5). The func- 
tional DNA markers for IcyE, though not consistent, were 
significantly associated only with lutein, and the ratio of 
(3 to a branch carotenoids, explaining 15 to 21% of the 
variations. However, IcyE-S TE did not significantly affect 
(3 to a branch ratio in the first year and lcyE-3 'indel was 
not associated with any of the traits. None of the markers 
for each gene had significant effects on total carotenoid 
content in both years which is consistent with previous 
association analyses results [12,34]. However, due to vari- 
ation in results of association and segregation mapping, 
earlier studies detected significant reduction of total carot- 
enoids for genotypes with favourable alleles of IcyE and 
crtRBl using segregating populations [34,35]. 



Combinatorial effects of lycE and crtRBl functional 
markers on carotenoid levels and profiles 

Fourteen unique genotypes were observed for IcyE and 
eight unique genotypes for crtRBl (Table 8). The two-way 
ANOVA combining the IcyE and crtRBl alleles revealed 
highly significant interaction effects for each carotenoid 
type and the derived traits. The combined effects of the 
alleles of the two genes were stronger than their separ- 
ate effects. The two genes model explained 38 to 89% of 
the total variation in carotenoid concentration. Individ- 
ual effects of the alleles were also highly significant for 
almost all carotenoids. The combined IcyE markers 
explained the least variation in the p-branch carotenoids 
(p-cryptoxanthine, p-carotene, zeaxanthine), while 
thecrtRBl markers explained the least variation in the 
a-branch carotenoid (lutein). The combined crtRBl 
markers had larger effects on individual and total pro- 
vitamin A carotenoids in comparison to the effects of 
the IcyE markers. 

Analysis of combinations of all of the six markers identi- 
fied 34 unique genotypes (Table 9). The vast majority of 
these genotypes were represented by only one inbred line 
each. The most common genotype, '0', G, 8|-, 1, 0, - |3 (cor- 
responding to /cj/£-5TE, -SNP (216), -3TnDel, crtRBl-5' 
TE, -indel4, 3TE; where the symbol '|' separates the two 
alleles of heterozygous loci, while the symbol represents 
any of the alternative alleles for the particular locus) was 
present in 49% of the inbred lines. The average estimated 
effect of the most frequent genotypes on beta-carotene 
was 3.54 \ig g" 1 , (which was 2.45 \xg g 1 less than the aver- 
age effect of all the genotypes and 6.9 \ig g" 1 less than the 
average effect of those genotypes containing the favourable 
alleles of both crtRBl 5TE and 3TE) (Table 9). Genotype, 
'0', G, 8|0, 2, 12, 1 contained the most optimal allelic 
combinations as it carried favourable alleles for 5 of the 6 
loci and was present only in one inbred line (Entry 107) 
derived from a backcross KU1409/SC55/KU1409 (Table 1, 
Table 4). Its estimated effects were 9.05 \ig g 1 for (3- 
carotene, 0.33 \ig g 1 for p-cryptoxanthine, 0.96 \ig g 1 for 
zeaxanthine and 13.29 \ig g l for total carotenoid. Al- 
though this genotype was predicted to be the best in terms 
of its allelic composition for the 6 markers, seven other 



Table 6 Linkage disequilibrium between markers of crtRBl and IcyE within the 130 yellow maize inbred lines 


R 2 * 


IcyE 5TE 


IcyE SNP (216) 


IcyE 3'indel 


crtRBl 5TE 


crtRBl indel4 


crtRBl 3TE 


IcyE 5'TE 




7.29E-1 1 


7.68E-13 


1 .09E-03 


3.11E-01 


8.50E-03 


IcyE SNP (216) 


0.3810 




3.92E-07 


1.43E-05 


5.66E-01 


1 .82E-04 


IcyE 3'indel 


0.3990 


0.2050 




1.53E-01 


1.13E-01 


1.01 E-01 


crtRBl 5'TE 


0.1010 


0.1880 


0.0201 




1.55E-01 


7.88E-19 


crtRBl indel4 


0.0115 


0.0041 


0.0208 


0.0166 




2.48E-01 


crtRBl 3'TE 


0.0658 


0.1350 


0.0262 


0.7570 


0.0144 





*Upper triangle contains p-value and lower triangle contains R 2 . 
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KU1409/SC55//KU 1409-98 (8.0)* 
KU1409/DE3//KU 1409-38 (14.6) 
KU1409/DE3//KU 1409-45 (15.4) 
KU1 409/DE3//KU1 409-49 (7.1 ) 
KU1409/DE3//KU 1409-53 (15.1) 
KU1409/DE3//KU 1409-43 (9.5) 
KU1409/DE3//KU 1409-55 (13.3) 
KU1409/DE3//KU1409-51 (8.4) 
KU1409/DE3//KU 1409-56 (9.4)* 
KU1409/DE3//KU 1409-44 (6.2) 
KU1409/DE3//KU 1409-48 (10.6) 
KU1409/DE3//KU 1409-99 (12.1) 
KU1409/DE3//KU 1409-39 (17.3) 
KU1409/DE3//KU 1409-47 (10.8) 
KU1409//KU1414-SR/A619-25 (9.4)* 
KU1409//KU1414-SR/A619-23 (11.6) 
KU1409//KU1414-SR/A619-24 (10.7) 
KU1409/NC358//KU1409-104 (7.9)* 
KU1409/SC55/KU1409-107 (9.6)* 
KU1409/NC358//KU1409-105 (4.0)* 
KU1409/NC358//KU1409-106 (13.0) 
SC55/KU1 41 4-SR//KU 1 41 4-SR-1 01 (1 0.9)* 
9450/KI21 -/DE3//9450/KI21 -1 1 4 (7. 1 ) 
KU 1 41 4-SR/CI7//KU 141 4-SR-1 23 (1 .7) 
9450/KI21 -/DE3//9450/KI21 -115 (9.8) 
9450/KI21 -/DE3//9450/KI21 -1 24 (7.9) 



Figure 2 Dendrogram of 26 inbred lines that have the best favourable alleles of crtRB1-5TE and crtRB-3TE marker. The pedigrees refer 
to the sources from which the inbred lines were derived. The numbers after the pedigrees are inbred line entry numbers. Numbers in parenthesis 
are mean (3-carotene concentration in ug" 1 DW. Entry 99 is the line used in [35] for developing segregating populations. Twenty of the inbred 
lines contained the favourable alleles of both markers except for those marked with * 



genotypes were found to be superior to this genotype in 
their estimated levels of p-carotene. The presence of an 
unfavourable IcyE insertion in the homozygous state did 
not alter the effect of this genotype significantly, based on 
the observation of the effects on carotenoids observed in 
another genotype '0', G, 8, 2, 12, 1 (N = 4), which lacked 
the lcyE-3' insertion. Genotype '0', G, 8, 2, 0, 1|- (N = 3) 
showed significantly better effect (p < 0.01) than the geno- 
type with the predicted best allelic composition ('0', G, 8|0, 
2, 12, 1), and had the strongest positive effect with an 
estimated average concentration of 15.03 \ig g" 1 for p- 
carotene and 15.08 \ig g" 1 for provitamin A. The major 
difference between the two genotypes is the lack of the 
favourable 12 bp insertion at crtRBl -indel4 in the 
former, which shows the negligible effect of this marker 
as was observed in the association analysis. Three inbred 
lines derived from KU1409/DE3/KU1409 contained the 
genotype '0', G, 8, 2, 0, 1 1 -. 

The favourable alleles of 5' and 3TE markers of crtRBl 
(alleles 2 and 1) were present in almost all genotypic com- 
binations that had large positive effects on p-carotene 
concentration ranging from 6.0 to 15.28 \ig g" 1 . The 



genotype 3, T, 8|0, 1, 0, 3 did not have any of the 
favourable alleles except the deletion allele representing 
/cy£-3'TE. Only one inbred line (Entry number 50) de- 
rived from a backcross KU1409/DE3/KU1409 carrying 
this genotype had estimated average effects of 7.56 \ig g' 1 
for p-carotene, 8.55 \ig g 1 for beta-cryptoxanthine, 
12.09 \ig g 1 for provitamin A and 43.3 (igg 1 for total ca- 
rotenoid. The total carotenoid concentration of this geno- 
type exceeded that of the average total carotenoid of 
those genotypes carrying the favourable alleles of crtRBl 
5'TE and 3'TE by 23.48 [ig g" 1 . This genotype also had 
7.32 \ig g" 1 higher p-cryptoxanthine than the average of 
those carrying the above mentioned allelic classes. These 
results were corroborated by the low ratio values of beta- 
carotene to p-cryptoxanthine, p-carotene to zeaxanthine 
and p-carotene to all carotenoids. Another exceptional 
genotype 3, G, 8, 1|2, 0, 3 containing the favourable 
allele of crtRBl 5TE in the heterozygous state showed 
a very weak effect on beta-carotene (1.05 \ig g" 1 ) and pro- 
vitamin A (1.65 \ig g 1 ) content. The weakest effect was 
detected from genotype '0,' G, 0, 1, 12, 3 (N = 3), which was 
devoid of the two best favourable alleles of crtRBl 5TE 
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Table 7 Marker-trait association of crtRBI and IcyE with carotenoid content of 130 yellow maize inbred lines 



Functional DNA (PCR) markers* 



Season 



Year 1 



Carotenoids* 



crtRBI 5TE 



crtRBI Indel4 



crtRBI 3TE 



IcyE 5TE 



IcyE 21 6(SNP) 



IcyE 3'lndel 



Year 2 



Combined* 





P 


R 2 P 


R 2 


P 


R 2 


P 


R 2 


P 


R 2 


Pear 


6.6E-13 


0.53 




1.2E-12 


0.51 










Pcry 


2.8E-08 


0.32 




1.5E-07 


0.28 










zeax 


3.7E-08 


0.31 




1.5E-07 


0.28 










lut 












3.1E-03 


0.11 


1 .6E-03 


0.08 


pva 


1 .2E-07 


0.26 2.3E-03 


0.09 


6.5E-07 


0.23 










InPcar/pcry 


8.6E-12 


0.47 




7.7E-1 1 


0.42 










InPcar/Zeax 


4.1E-11 


0.45 




1.7E-10 


0.42 










InPcar/all 


5.4E-08 


0.30 




3.7E-08 


0.30 










Inpbr/abr 


1.3E-04 


0.15 




4.7E-04 


0.13 






6.7E-04 


0.10 


Pear 


4.4E-10 


0.38 




7.9E-1 1 


0.42 










Pcry 


4.6E-08 


0.30 




7.7E-08 


0.29 










zeax 


1 .2E-05 


0.20 




9.5E-06 


0.20 










lut 












1 .2E-04 


0.18 


2.7E-05 


0.15 


tpva 


2.2E-05 


0.17 




4.8E-06 


0.20 










InPcar/pcry 


2.1E-14 


0.63 




1.3E-13 


0.58 










InPcar/Zeax 


8.1E-11 


0.44 




1.3E-11 


0.47 










InPcar/all 


5.4E-09 


0.34 




2.5E-10 


0.41 










Inpbr/abr 












2.6E-04 


0.16 


9.1E-07 


0.21 


Pear 


3.4E-12 


0.49 




2.5E-12 


0.50 










Pcry 


1.1E-08 


0.34 




3.2E-08 


0.31 










zeax 


1 .6E-06 


0.24 




1 .8E-06 


0.23 










lut 












4.3E-05 


0.20 


1.5E-06 


0.20 


tpva 


1 .2E-08 


0.31 




6.3E-09 


0.33 










InPcar/pcry 


1.8E-14 


0.63 




8.4E-14 


0.59 










InPcar/Zeax 


1.9E-11 


0.47 




1.1E-11 


0.48 










InPcar/all 


9.3E-09 


0.33 




4.5E-10 


0.39 










Inpbr/abr 


8.1E-03 


0.08 








7.2E-04 


0.14 


2.0E-06 


0.20 



6.9E-04 0.12 



5.8E-05 0.17 



8.3E-03 0.08 



^Carotenoid abbreviations are described under Table 3. pbr/abr = p-branch/a-branch carotenoids, 

pcar/pcry = p-carotene/p-cryptoxanthine, pcar/Zeax = p-ca rote ne/zeaxan thine; pear/all = p-carotene/total carotenoids. 

** Only significant values are indicated. Names of functional markers for each gene are described in Table 2. Significance thresholds are 1.67E-04 and 8.3E-3 at 1% 
and 5%, respectively, after Bonferroni multiple test correction (alpha/number of markers analyzed per trait, i. e., 0.01/6= 1.7E-04 and 0.05/6 = 8.3E-3) (TASSEL user's 
guide 2011). 

***BLUEs calculated from the two year data illustrated using bar graph in Figure 3. 



and 3TE (Table 4), which had relatively low level of total 
carotenoids (10.18 Lig g" 1 ). Overall, the average effects of 
the genotypes harbouring both favourable alleles of crtRBI - 
5'TE and -3TE (N = 23) resulted in 7.2 Lig g" 1 increase or 
3.22 fold increase in (3-carotene as compared to the effects 
of genotypes without any of the favourable alleles (N = 
103). The reduction in total carotenoid between the two 
sets of allelic classes was found to be negligible (from 23.5 
to 18 Lig g 1 ). 

Discussion 

Maize is an important staple food crop for food and liveli- 
hood security in Africa. Since the 1970s, IITA has had a 



maize breeding program to develop tropical maize lines 
that are high-yielding and adapted to growing conditions 
across Africa. To help alleviate micronutrient deficiencies 
amongst the poor whose diets are highly dependent on 
maize, the development of tropical adapted maize lines 
with elevated levels of carotenoids, in particular provitamin 
A, is a major maize improvement goal. 

To facilitate the development of high vitamin A trop- 
ical maize varieties, the carotenoid profile and content 
of 130 tropical adapted inbred lines within the conven- 
tional breeding program for tropical maize at IITA were 
analysed over two seasons under field conditions. Within 
these 130 inbred lines, 28 lines were found to have a 
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□ Pear 

■ fkry 

■ zeax 

■ lut 

□ pva 

■ ln(3car/(3cry 

□ ln(3car/Zeax 

■ ln(3car/all 

■ ln(3br/abr 



crtRBl-5'TE crtRBl 3'TE lcyE-5TE lcyE-216 (SNP) lcyE-3'lndel 
Name of functional markers 

Figure 3 Significant marker-trait associations estimated using BLUEs calculated based on the two year field data. Significance thresholds: -log 
(1 .7E-04) = 4 at a = 0.01 and -log (8.3E-3) = 2 at a = 0.05. Carotenoid names abbreviations described in Tables 3 and 7. R 2 Values indicated on top of each 
bar represent percent variation explained. 



provitamin A content of 8.0 to 17.25 \ig g" 1 which is higher 
than previously reported by [13] and is comparable to the 
highest provitamin A level described in [14]. This finding 
highlights the importance of introgressing the best 
favourable alleles of provitamin A from temperate germ- 
plasm into tropical adapted inbred lines. The maximum 
average estimated level of total carotenoids detected was 
42.71 ug g" 1 , which was much lower than the 100 ug g" 1 
previously reported [11]. The identification of inbred lines 
with high total carotenoid levels is considered to be useful 
if the high influx of substrates to the carotenoid biosyn- 
thesis pathway favors an increase in provitamin A caroten- 
oids in maize endosperm [32]. 

DNA markers detecting polymorphisms in genes that 
are functionally responsible for changes in phenotypes 
can be called functional markers [27,45]. Eight functional 
markers of three key carotenoid genes PSY1 [12], IcyE [33] 
and crtRBl [34], (previously developed using different as- 
sociation panels of diverse temperate, sub-tropical and/or 
tropical yellow maize inbred lines), were considered for 
validation in this study. The functional markers for the 
two genes IcyE and crtRBl have been described as one of 
the most exciting discoveries for maize endosperm pro- 
vitamin A improvement endeavors [14]. However, inde- 
pendent studies demonstrated some inconsistencies in the 
effects of these markers [11,35,36] which necessitates 



additional investigations so that such markers can be de- 
ployed in maize breeding in a more robust and predictable 
manner. In the earlier studies, the panels and populations 
used for developing and validating the functional markers 
were largely of temperate origin and had low frequencies 
of the favourable alleles of the most significant markers. 
In our field study over two years, these functional markers 
have been analysed for their efficacy in diverse maize 
inbred lines derived mainly from populations containing a 
mixture of tropical and temperate germplasm in their 
pedigrees. Our study clearly demonstrates that the effects 
are heritable thereby can facilitate the development of 
robust maize varieties with elevated provitamin A levels. 

The functional markers for PSY1 were monomorphic 
for the favourable allelic variants across all the inbred 
lines possibly because of the highly conserved nature of 
the PSY1 gene within and across species [46]. For in- 
stance, Fu et al. [33] have observed that the favourable 
alleles of PSY1 were fixed within the tropical genetic 
background of the panels used in their study. The vari- 
ation in total carotenoid content observed in our study 
could be due either to the presence of some rare func- 
tional variation within the PSY1 gene and its regulatory 
regions (of a genetic or epigenetic nature) and/or other 
"modifier" genes that are involved in the carotenoid 
biosynthesis in the genotypes tested [32]. 



Table 8 Phenotypic variation explained (R 2 ) by individual and combined effects of IcyE and crtRBl alleles 
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3 
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0.42 
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3 


10 


0.24 


0.35 


0.28 


0.63 


0.48 


NS 


0.23 


0.71 


0.60 


0.50 



*Names and symbols of polymorphic sites and their allelic variants are described in Table 2. **R from combined ANOVA computed for each of the 3 haplotype 
scenarios across 2 years. NS, non-significant at alpha 0.5, the rest are significant at alpha < 0.0001. Abbreviation of carotenoids' names described under Table 2 and 3. 
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-0.67 ± 0.36 


-1.68 ±0.52 
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N = number of inbred lines, NA = Not available (missing marker score). SE = Standard error of least square means. Symbol '\ separates heterozygous alleles of a locus. SNP (216) is a dominant marker for the T allele, 
thus T can be 'TP or TG'. 
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The markers for IcyE and crtRBl were polymorphic 
across the maize inbred lines, where in the 3' and 5'TE 
markers of crtRBl exhibited strong association with vari- 
ation in p-carotene content of the inbred lines. In contrast, 
the effects of IcyE markers were found to be weak and in- 
consistent in the present study, which was in line with pre- 
vious results [11,35]. The germplasm within the tropical 
maize gene pool is known to be more diverse than that in 
the temperate maize gene pool. In a previous study Yan 
et al. [34] detected the favourable allele of crtRBl -5'TE 
only in the temperate yellow maize germplasm, with a fre- 
quency of less than 3%. In our study on tropical maize 
germplasm, this allele occurred at a relatively high fre- 
quency of 18%. The high linkage disequilibrium (R 2 = 0.76) 
between the 3TE and 5TE polymorphisms of crtRBl de- 
tected in our study deviates from a previous report [34] 
that found no linkage disequilibrium (R 2 = 0.02). It is prob- 
able that the two favourable alleles were introgressed (from 
temperate donor inbred lines into the tropical adapted ma- 
terials) together as genetically linked alleles that led to the 
observed strong linkage disequilibrium between the two 
markers. Such linkage disequilibrium makes the estimation 
of the independent effect of each marker difficult. How- 
ever, based on results of previous studies [34,35], it can be 
argued that both markers could be contributing to the 
strong association of crtRBl with provitamin A in maize 
endosperm, with the 5TE polymorphism contributing the 
largest effect. 

One of the maize lines (derived from a back cross 
involving the maize line DE as a donor parent) was also 
used for developing five different segregating populations 
used by [35] to test the effect of crtRBl-3'TE polymorph- 
ism [35]. In our study, this line along with 25 other inbred 
lines carried the best favourable alleles of crtRBl -5TE 
and -3TE polymorphisms and was used to evaluate to 
what extent the lines carrying the two favourable alleles in 
this study were different from the segregating population 
used in [35]. UPGMA cluster analysis of SNPs based on 
Neis 1972 distance (Figure 2) separated the lines accord- 
ing to their genetic backgrounds. Substantial genetic vari- 
ation was found among the lines originating from the 
same backcross. The line used by [35] was clustered with 
one of the major groups thus underpinning the diversity 
of the lines used in our study. In previous analyses, the ef- 
fect of crtRBl-5'TE, the marker that was reported to have 
the largest effect in the work of Yan et al. [34], was not 
reported in the validation study of [35]. Hence, our study 
fills a major gap by now providing the marker-trait associ- 
ation results for crtRBl -5'TE in tropical adapted maize 
germplasm. 

Almost all of the inbred maize lines with high levels of 
p-carotene and total provitamin A carried the favourable 
alleles of the most significant functional markers of 
crtRBl-3 Y and crtRBl-5' TE. An inconsistency detected 



in our marker-trait study was a maize inbred line that 
showed unexpectedly low p-carotene and provitamin A 
(<2.0 \xg g" 1 DW) although it carried the favourable allele 
of crtRBl -5TE polymorphism. This inbred line also had 
relatively lower total carotenoid content (14.14 \ig g" 1 
DW), possibly suggesting that high influx of substrates into 
the carotenoid biosynthesis pathway may be an important 
factor to realize the desired action of the favourable alleles 
[11,33,34]. Hence, introgression of these favourable alleles 
into adapted maize germplasm with high total carotenoid 
content could be a strategy for increasing levels of p- 
carotene and total provitamin A. 

The combinatorial analysis of the functional polymor- 
phisms for the two genes revealed larger effects than 
those observed for alleles of each gene independently. 
This finding is in agreement with results of previous 
studies [34,35]. In particular, our analysis identified a 
number of superior genotypes of maize inbred lines that 
have carotenoid levels of relevance to provitamin A level 
enhancement that exceed previously obtained levels 
within IITAs breeding program. Our results indicates 
that the functional markers for crtRBl markers have the 
strongest potential to accelerate genetic gain for en- 
hanced p-carotene content in tropical maize breeding 
programs [35]. Given that the carotenoid biosynthesis 
pathway is most likely conditioned by a number of genes 
and regulatory elements, and that all of the variation in 
provitamin A levels are not accounted for, it can be 
worthwhile to consider genomic selection approaches for 
enhanced provitamin A carotenoid levels. 

Conclusions 

The first generation of provitamin A enriched maize hy- 
brids have been developed for Nigeria and Zambia re- 
cently, which will most likely spread to other African 
countries with similar agro-ecologies over the coming years 
[47,48]. In this two year field study on tropical maize 
germplasm we have demonstrated the strong associa- 
tion between favourable alleles detected by the crtRBl- 
5'TE and -3'TE functional markers and high levels of 
P-carotene. Our study found these two markers to be in 
strong linkage disequilibrium in the tropical maize germ- 
plasm, raising the possibility that one of the polymorphic 
sites (e.g. the 5'TE marker) could be targeted to reduce 
costs associated with PCR genotyping. The high provitamin 
A inbred lines harbouring combinations of the favourable 
alleles of the crtRBl -5'TE and crtRBl -3'TE markers can be 
used to speed up the development of the next generation 
of high provitamin A tropical maize hybrids for production 
in Sub-Saharan Africa, thus contributing to the alleviation 
of hidden hunger due to vitamin A deficiency. 

Competing interests 

The authors declare that they have no competing interests. 



Azmach et al. BMC Plant Biology 2013, 13:227 
http://www.biomedcentral.com/1471-2229/13/227 



Page 15 of 16 



Authors' contributions 

GA conducted the field and laboratory experiments, carried out the statistical 
analysis, and prepared the draft of the manuscript. AM contributed to the 
conception of the project, supplied maize inbred lines developed by AM, 
supervised the field work and revised the manuscript. MG closely supervised 
the laboratory work in IITA and revised the manuscript. CS contributed to 
the conception of the project, securing of funding, participated in the 
supervision of the overall work and contributed to the drafting and 
finalization of the manuscript. CS, AM and MG participated in the PhD 
supervision of GA. All authors read and approved the final manuscript. 



Acknowledgements 

This work was part of a collaborative PhD project between the International 
Institute for Tropical Agriculture (IITA) and the Plant and AgriBiosciences 
Research Centre (PABC) at the National University of Ireland Galway (NUI 
Galway). The research was financed by Irish Aid, through IITA and NUI 
Galway. We thank the staff of Maize Improvement Unit and Bioscience 
Centre of IITA for their support in conducting the experiments. 

Received: 20 September 2013 Accepted: 19 December 2013 
Published: 28 December 2013 



References 

1. Cunningham FX, Gantt E: Genes and enzymes of carotenoid biosynthesis 

in plants. Annu Rev Plant Physiol Plant Mol Biol 1998, 49:557-583. 

2. Fraser PD, Bramley PM: The biosynthesis and nutritional uses of 

carotenoids. Prog Lipid Res 2004, 43:228-265. 

3. Krinsky Nl, Johnson EJ: Carotenoid actions and their relation to health 
and disease. Mol Aspects Med 2005, 26:459-516. 

4. West KP Jr, Darnton-Hill I: Vitamin a deficiency. In Nutr heal dev ctries SE - 13. 
Edited by Semba RD, Bloem MW.Totowa, NJ: Humana Press; 2008:377-434. 

5. Cuttriss AJ, Cazzonelli CI, Wurtzel ET, Pogson BJ: Carotenoids. In Adv hot res 
biosynth vitam plants part A vitam A, Bl, B2, B3, B5, Volume 58. Edited by 
Rebeille F, Douce R. Amsterdam, The Netherlands: Elsevier Ltd. Academic 
Press; 2011:1-36. 

6. World Health Organisation: Vitamin a. In Vitam A miner requir hum nutr. 
Geneva, Switzerland: WHO; 2004. 

7. Sherwin JC, Reacher MH, Dean WH, Ngondi J: Epidemiology of vitamin a 
deficiency and xerophthalmia in at-risk populations. Trans R Soc Trop Med 
Hyg 20U, 106:205-14. 

8. WHO: Global prevalence of vitamin A deficiency in populations at risk 1995-2005. 
Geneva: World Health Organization; 2009. 

9. Darnton-Hill I, Webb P, Harvey PWJ, Hunt JM, Dalmiya N, Chopra M, Ball MJ, 
Bloem MW, de Benoist B: Micron utrient deficiencies and gender: social 
and economic costs. Am J Clin Nutr 2005, 81:11 98S-1 205S. 

10. The World Bank: Repositioning nutrition as central to development A strategy 
for large-scale action. Washington DC: The World Bank; 2006:246. 

1 1 . Burt AJ, Grainger CM, Smid MP, Shelp BJ, Lee EA: Allele mining of exotic 
maize germplasm to enhance macular carotenoids. Crop Sci 201 1, 
51:991-1004. 

12. Harjes CE, Rocheford TR, Bai L, Brutnell TP, Vallabhaneni R, Williams M, 
Wurtzel ET, Kandianis CB, Sowinski SG, Stapleton AE, Yan J, Buckler ES: 
Natural genetic variation in lycopene epsilon cyclase tapped for maize 
biofortification. Science 2008, 319:330-3. 

13. Menkir A, Liu W, White WS, Maziya-Dixon B, Rocheford T: Carotenoid 
diversity in tropical-adapted yellow maize inbred lines. Food Chem 2008, 
109:521-529. 

14. Pixley K, Rojas NN, Babu R, Mutale R, Surles R, Simpungwe E: 
Biofortification of maize with provitamin a carotenoids. In Carotenoids 
hum heal SE - 17. Edited by Tanumihardjo SA. LA - English: Humana Press; 
2013:271-292. 

15. Bouis HE, Hotz C, McClafferty B, Meenakshi JV, Pfeiffer WH: Biofortification: 
a new tool to reduce micronutrient malnutrition. Food Nutr Bull 201 1, 
32(1):S31-40. 

16. Wurtzel E, Cuttriss A, Vallabhaneni R: Maize provitamin a carotenoids, current 
resources, and future metabolic engineering challenges. Front Plant Sci 
2012, 3:1-12. 

17. Forum for Agricultural Research in Africa: Patterns of change in maize 
production in africa: implications for maize policy development. Accra, Ghana: 
FARA; 2009:1-8. 



18. Smale M, Byerlee D, Jayne T: Maize revolutions in sub-saharan africa. In An 
African green revolut SE - 8. Edited by Otsuka K, Larson DF. Amsterdam, 
Netherlands: Springer; 2013:165-195. 

19. Bouis HE, Welch RM: Biofortification-a sustainable agricultural strategy for 
reducing micronutrient malnutrition in the global south. Crop Sci 2010, 
50(1):S-20-S-32. 

20. Meenakshi JV, Johnson NL, Manyong VM, Degroote H, Gonzalez C, Garci J, 
Meng E: How cost-effective is biofortification in combating micronutrient 
malnutrition? An ex ante assessment. World Dev 2010, 38:64-75. 

21. Msangi S, SulserTB, Bouis A, Hawes D, Batka M: Integrated economic 
modeling of global and regional micronutrient security. Washington DC: 
Harvest Plus Working Paper no. 5, International Food Policy Research 
Institute (IFPRI); 2010:75. 

22. Qaim M, Stein AJ, Meenakshi JV: Economics of biofortification. Agric Econ 
2007, 37:119-133. 

23. Underwood BA: Overcoming micronutrient deficiencies in developing 
countries: is there a role for agriculture? Food Nutr Bull 2000, 
21:356-360. 

24. Pfeiffer WH, McClafferty B: Harvestplus: breeding crops for better 
nutrition. Crop Sci 2007, 47:S-88. 

25. Mishra P, Singh NK: Spectrophotometry and TLC based characterization 
of kernel carotenoids in short duration maize. Maydica 2010, 55:95-100. 

26. Moose SP, Mumm RH: Molecular plant breeding as the foundation for 
21st century crop improvement. Plant Physiol 2008, 147:969-977. 

27. Andersen JJR, Lubberstedt T: Functional markers in plants. Trends Plant Sci 
2003, 8:554-560. 

28. Cazzonelli CI, Pogson BJ: Source to sink: regulation of carotenoid 
biosynthesis in plants. Trends Plant Sci 2010, 15:266-274. 

29. DellaPenna D, Pogson BJ: Vitamin synthesis in plants: tocopherols and 
carotenoids. Annu Rev Plant Biol 2006, 57:71 1-738. 

30. Farre G, Sanahuja G, Naqvi S, Bai C, Capell T: Travel advice on the road to 
carotenoids in plants. Plant Sci 2010, 179:28-48. 

31. Hirschberg J: Carotenoid biosynthesis in flowering plants. Curr Opin Plant 
Biol 2001,4:210-218. 

32. Vallabhaneni R, Wurtzel E: Timing and biosynthetic potential for 
carotenoid accumulation in genetically diverse germplasm of maize. 
Plant Physiol 2009, 150:562-572. 

33. Fu Z, Chai Y, Zhou Y, Yang X: Natural variation in the sequence of psyl 
and frequency of favorable polymorphisms among tropical and 
temperate maize germplasm. Theor Appl Genet 2013, 126:923-935. 

34. Yan J, Kandianis C, Harjes C, Bai L: Rare genetic variation at zea mays crtrbl 
increases -carotene in maize grain. Nat Genet 2010, 42:322-327. 

35. Babu R, Rojas NP, Gao S, Yan J, Pixley K: Validation of the effects of 
molecular marker polymorphisms in Icye and crtrbl on provitamin a 
concentrations for 26 tropical maize populations. Theor Appl Genet 2013, 
126:389-99. 

36. Vignesh M, Hossain F, Nepolean T, Supradip S, Agrawal PK, Guleria SK, 
Prasnna BM, Gupta HS: Genetic variability for kernel (3 -carotene and 
utilization of crtrbl 3' te gene for biofortification in maize (zea mays I.). 
Indian J Genet Plant Breed 2012, 72:189-194. 

37. Howe JA, Tanumihardjo SA: Carotenoid-biofortified maize maintains 
adequate vitamin a status in mongolian gerbils. J Nutr 2006, 
136:2562-2567. 

38. US Institute of Medicine: Dietary reference intakes for vitamin A, vitamin K, 
arsenic, boron, chromium, copper, iodine, iron, manganese, molybdenum, 
nickel, silicon, vanadium, and zinc. Washington, DC: The National Academies 
Press; 2001. 

39. Saxton A: A macro for converting mean separation output to letter 
groupings in proc mixed. Proc 23rd SAS Users Gr Int 1998:1246-1246. 

40. Yu J, Pressoir G, Briggs W, Bi I, Yamasaki M: A unified mixed-model method 
for association mapping that accounts for multiple levels of relatedness. 
Nat Genet 2006, 38:203-208. 

41. Bradbury PJ, Zhang Z, Kroon DE, Casstevens TM, Ramdoss Y, Buckler ES: 
Tassel: software for association mapping of complex traits in diverse 
samples. Bioinformatics 2007, 23:2633-2635. 

42. Elshire R, Glaubitz J, Sun Q: A robust, simple genotyping-by-sequencing 
(gbs) approach for high diversity species. PLoS One 201 1, 6:1-10. 

43. Liu K, Muse SV: Powermarker: an integrated analysis environment for 
genetic marker analysis. Bioinformatics 2005, 21:2128-9. 

44. Nei M: Genetic distance between populations. Am Nat 1972, 
106:283-292. 



Azmach et al. BMC Plant Biology 2013, 13:227 
http://www.biomedcentral.com/1471-2229/13/227 



Page 16 of 16 



45. Varshney RK, Graner A, Sorrells ME: Genomics-assisted breeding for crop 
improvement. Trends Plant Sci 2005, 10:621-630. 

46. Palaisa K, Morgante M: Contrasting effects of selection on sequence 
diversity and linkage disequilibrium at two phytoene synthase loci. 

Plant Cell 2003, 15:1795-1806. 

47. IITA. http://www.iita.Ofg. 

48. Harvestplus. http://www.haFvestplus.org. 



doi:1 0.1 186/1471-2229-13-227 

Cite this article as: Azmach et al.: Marker-trait association analysis of 
functional gene markers for provitamin A levels across diverse tropical 
yellow maize inbred lines. BMC Plant Biology 2013 13:227. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at (^\ RiftMM i rpntral 

www.biomedcentral.com/submit \^ ™omea centra I 



